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Abstract 

We study collider signatures for extra scalar bosons in the inert doublet model at the interna- 
tional linear collider (ILC). The inert doublet model is a simple extension of the standard model 
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stable and a candidate of the dark matter as long as it is electrically neutral. Taking into account 
the constraints from various theoretical and phenomenological conditions, we perform a simulation 
study for the distinctive signatures of the extra scalars over the standard-model background con- 
tributions at the ILC with the center-of-mass energy of \fs = 250 GeV and 500 GeV. We further 
discuss observables for determination of the mass of the scalars. We find that the parameter regions 
which cannot be detected at the large hadron collider can be probed at the ILC. 
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I. INTRODUCTION 



In July 2012, a Higgs-like particle was found at the large hadron collider (LHC) JT|, 0. 
The particle looks like mostly a Higgs boson in the standard model (SM), but the detail 
properties of the particle and the whole structure of the Higgs sector have not yet been 
revealed. It is widely believed that the SM has to be extended, since it cannot explain 
the dark matter, (tiny) neutrino masses, and the baryon asymmetry in the Universe, etc. 
Although in the Higgs sector of the SM, only one SU(2)^-doublet scalar field is introduced, 
there is no theoretical guideline for this choice. Thus, the Higgs sector may be a solid target 
to probe new physics beyond the SM. 

The inert doublet model (IDM) is one of the simplest extensions of the SM, where an 
additional SU(2)i-doublet scalar field is introduced, which is odd under the unbroken Z2 
symmetry [|3], |J. As in the case in the general two Higgs doublet model, four kinds of addi- 
tional scalars appear as physical states, namely neutral CP-even state (if), neutral CP-odd 
state (A) and charged scalar states (H ), all of which are called inert scalars. Due to the Z2 
symmetry, Yukawa interactions of the inert scalars to SM fermions are forbidden, and the 
possible flavor-changing neutral current is absent at the tree level. The lightest inert particle 
(LIP) is stable, because of the ^-parity conservation. Therefore, the model provides a dark 
matter candidate In addition to that, the model has rich phenomenological features 



such as the electroweak symmetry breaking |T]J , electroweak phase transition ||i2HT5|l , radia- 



tively generating neutrino masses by introducing Z 2 -odd right-handed neutrinos [16], and 



leptogenesis |L7|, pl|, etc. 

Collider signatures of the inert scalars in the IDM have been studied in the litera- 
ture [|], ||, [10], p!^-p2|. In Ref. PC|| , bounds on the inert scalar masses are obtained by 



using the experimental results at the LEP II [p^-p5|. Since the inert scalars do not have 
QCD interactions, it is not suited to search for them at hadron colliders. Even though 
the parameter regions where the inert scalars could be discovered at the LHC are pointed 
out H 21, 22], detailed analysis on these scalars such as the precise determination of these 
masses and quantum numbers would be difficult. 

In this letter, we study collider phenomenology for the inert scalars at the international 
linear collider (ILC). As it is a machinery for precision measurements of Higgs boson proper- 
ties, the extended Higgs sector can also be investigated in details. We study the characteristic 



signatures, corresponding backgrounds and the kinematical observables in the processes of 
HA associated production as well as H + H~ pair production. Earlier studies can be found 
e.g. in Refs. |26], |27j where the charged scalar pair production process is studied in detail, 
and in Ref. [2S| where HA associated production is also studied. On the other hand, our 
study includes all the available processes and decay modes, and the simulation analysis for 
the signal and the background contributions with appropriate kinematical cuts. Further- 
more, we also discuss the procedure for the mass determination of the inert scalars which 
can be performed at an early stage of the experiment, and introduce a new observable which 
could be useful to determine the inert scalar masses more precisely than the variables used 



in Refs. 27, 28 



The rest of the letter is organized as follows. In Sec. [II], we review the inert doublet 
model and introduce the benchmark points used in our simulation study. Then we present 
the simulation studies for observing the inert scalars and determining the masses of them 
in Sec. [TTT|. Sec. [TV] and Sec. [V] are devoted to discussions and conclusions, respectively. In 
Appendix, we evaluate a new observable at e + e~ colliders which is used in our analysis for 
the mass determination. 



II. THE INERT DOUBLET MODEL 

In the IDM, the two scalar doublet fields, $1 and $2, have even and odd parity under 
the Z 2 symmetry, respectively. The most general scalar potential can be written as 

V{$ u $ 2 ) = fil l^l 2 + ll\ |$ 2 | 2 + y |$i| 4 + y |$ 2 | 4 + A 3 |$i| 2 |$ 2 | 2 



+ \ A 



2 + {^ 2 



+ H.c. , (1) 



with seven real parameters (/if, Ai, A2, A3, A4, A5). We note that the potential is invariant 
under the CP transformation. The potential has to satisfy theoretical constraints, such 
as the vacuum stability || |{| and the perturbativity Q. By the vacuum stability at the 
tree level, the quartic terms are constrained as Ai > 0, A 2 > 0, \J\\\<i + A 3 > 0, and 



a/AiA 2 + A3 + A4 — I A5 1 > 0. We consider the case where n\ < 0, Ai/i| > A3//1 and 
Ai/x| > (A3 + A4 + |A 5 |)/if are satisfied ||, so that <3> 2 does not acquire the vacuum expectation 



value (VEV) and only $1 plays a role of the "Higgs-boson" . By denoting 

(2) 





where v is the VEV, v = a/— 2/xf/Ai(~ 246 GeV), the masses of these scalars are expressed 
as 

ml = \iv 2 , (3a) 
m H+ = l4 + ^v 2 , (3b) 
m \ — A*2 + ^(^3 + A 4 + A 5 )t> 2 , (3c) 
™a = A*l + 7^3 + A 4 - X 5 )v 2 . (3d) 

Thus, the seven parameters in the Higgs potential can be replaced by the VEV v, four masses 
of the Higgs boson and inert scalars, (m h ,mH+,mH, m A), the scalar self-coupling constant 
A2, and Aff(= A3 + A4 + As) for example. To force the LIP to be electrically neutral, so that 
it can be a candidate of the dark matter, A4 < \X^\ must be satisfied f29fl . Depending on the 
sign of A5, either H or A becomes the LIP. Since phenomenological constraints and collider 
signatures are exchangeable between the two cases, we take H as the LIP (A5 < 0) hereafter. 

In this letter, we consider four benchmark points for the masses of inert scalars listed in 
Table |, which satisfy all the available theoretical and also phenomenological constraints. 
Phenomenological constraints known so far |3(J consist of the direct searches at the LEP II 
experiments |2(J , electroweak precision tests j| , and the relic abundance and direct detection 



of the dark matter [§,0. In the four benchmark points, the mass of H is fixed to 65 GeV, 
so that it does not induce the invisible decay of the SM Higgs boson. While it could be 
up to ~ 80 GeV concerning the dark matter relic abundance analysis |5]-|TU|], the collider 



phenomenology does not change qualitatively by varying it in that region. Above that value, 
HH — > W + W~ process gives a too large annihilation cross section, so that the predicted 
relic density in the IDM results into a lower value than that from the WMAP experiment. 
Below 80 GeV, the processes HH — > ff can give proper annihilation cross sections using the 
Higg boson mediated diagram, at the same time with avoiding the constraints from direct 
searches by choosing the hHH coupling constant A# as an appropriate value. Here, we 
have assumed the mass of the Higgs boson to m h = 126 GeV. We study the case where the 



masses of H and A are close to each other (I, III) for which the LEP and LHC experiments 



can not probe || |19H2"2|j . The other cases are when uia — van is medium (II) or large such 
that the Z-boson from A — > HZ becomes on-shell (IV). For the VT-bosons in H — > W ± H, 
we consider the off-shell (I, II) and on-shell (III, IV) cases. These mass differences result 
into additional contributions [[| to the electroweak S and T parameters which are consistent 
with the current experimental data BT . 



A2 and Xh do not enter in our collider analysis. Xh can be determined by the relic 
abundance which should be consistent with the recent WMAP search result || and also 
the bounds on the cross section of dark matter direct production. On the other hand, the 
method to determine A2 has not yet been established. 

For our four benchmark points, the production cross sections of inert scalars at the ILC 
are large enough to be tested. In Table d we also list the cross sections of HA production 
and H + H production at y/s = 250 GeV and 500 GeV. The cross section of HA production 
takes the largest value, i.e. 186 fb at yfs = 190 GeV, 78 fb at yfs = 280 GeV, and 46 fb at 
y/s = 350 GeV in the cases (I, III), (II), and (IV), respectively. The cross section of H + H~ 
production takes the largest value, i.e. 96 fb at a/s = 380 GeV and 53 fb at a/s = 500 GeV 
for the cases (I, II) and (III, IV), respectively. At y/s = 1 TeV, the production cross sections 
are about 10 fb and 20 fb for HA production and H + H~ production, respectively, for all 
the four benchmark points. For the cases (II, IV), H decays into W ± H predominantly, 
where we admit the VT-boson to be off-shell if m#± — niu < rriw- While for the cases (I) 
and (HI), H — > W ± A decay would be sizable as well, with the branching ratios about 32% 
and 27%, respectively. The decay of the A-boson is dominated by A — > Z^H. 

There are two important effects in the property of the SM-like Higgs boson, the invisible 





Inert scalar masses 
m H [GeV] m A [GeV] m H ± [GeV] 


ILC cross sections [a/s = 
a e + e -^HA [fb] <?V< 


250 GeV (500 GeV)] 
e-^H+H- [fb] 


(I) 


65. 


73. 


120. 


152. (47.) 


11. (79.) 


(II) 


65. 


120. 


120. 


74. (41.) 


11. (79.) 


(III) 


65. 


73. 


160. 


152. (47.) 


0. (53.) 


(IV) 


65. 


160. 


160. 


17. (35.) 


0. (53.) 



TABLE I. Masses of inert scalars and ILC cross sections for our four benchmark points. 
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decay [19|, Q and the charged scalar loop contribution to the two-photon decay ampli- 
tude JED), [33], [34]]. The invisible decay mode opens if > Iran- If this is the case, the 
branching ratio for the invisible decay mode can be typically several tens percent [ |T9| [33 . 
Thus, it could be discovered at the LHC [Q. The two-photon branching ratio of the SM-like 
Higgs boson can be directly related to _R 77 = a(pp — > h — > 77)/c(pp — > h — > 77)sm, because 
the production cross-section of the Higgs-boson at the LHC is not modified in the IDM. It is 
shown |33| that _R 77 can be enhanced with relatively light H ± (m#± < 130 GeV), negative 
A3 and no invisible decays. 



III. COLLIDER SIGNATURES IN THE IDM AT THE ILC 



In this section, we perform the simulation studies for the detection and the mass deter- 
mination of the inert scalars in the IDM at the ILC. The kinematical distributions are cal- 
culated by using MadGraph [[35] at the parton level with basic cuts for event generations. For 



charged leptons, we set p? > 1 GeV, \r)i\ < 2.5 and ARa > 0.2 for the isolation requirement, 
where px is the transverse momentum, r] is the pseudo-rapidity, and AR(= a/ Arf + A<p 2 ) 
is the distance of the two particle in the i] — <fi plane. For jets, we restrict ourselves with 
p 3 T > 5 GeV, |?7j| < 2.5 and ARjj, ARij > 0.4. Furthermore, we require the missing trans- 
verse energy Ifi? to be greater than 10 GeV to reduce background events from two photon 
scattering processes. The two photon scattering processes and QED radiation effects are 
not estimated in our analysis. 

We note that, these cuts may be rather conservative, so that our parton level analysis 
makes sense. At e + e~ colliders, hadronic final-states would be utilized even if they do not 
form narrow jets. Then, the cuts on pj- and the isolation for partons may not be necessary. 
We include these cuts so that the number of jets in an event can be more easily accounted 
the number of outgoing partons in the processes. In case we loose these cuts, the number 
of available events would increase, but we expect more background contributions from the 
events with less partons. An isolation requirements for leptons may be weakened as well in 
the real experiment. 



A. e + e — > HA process 

First, we consider the HA production process followed by A — > Z^H decay. Since H is 
neutral and stable, it escapes from detection. Thus, it gives the signature with a dilepton 
(dijet) plus large missing energy. Expected background contributions come from dilepton 
(dijet) plus two neutrinos production in the SM. 

First we study the collider signature of these events at \fs = 250 GeV, for the parameter 
sets (I, III) and (II) where A decays into H and off-shell Z-boson. To reduce the SM 
background contributions, we apply following kinematical cuts; the scaled acoplanarity aC op, 
which is the acoplanarityQ multiplied by the sine of the smallest angle between a lepton (jet) 
and the beam axis, is larger than 100°; | cos6*«| < 0.8 for dilepton or | cos 6^1 < 0.6 for dijet, 
where 9u(jj) is the polar angle of the dilepton (dijet) 3-momenta. 

In the top left (right) panel in Fig. [I], we see the dilepton (dijet) energy distributions for 
the signal and background processes at yfs = 250 GeV with the integrated luminosity of 
£ int = 250 fb _1 . After the cuts described above, background events are well reduced. The 
difference of the number of the signal events in dilepton and dijet signature comes from the 
branching ratio of the Z-boson, and that of the background events comes from the absence 
of WW production in the dijet case. 

The endpoints of the -EWjjj) distribution are related with the masses of H and A as 



UM 4 V m\ 



s s ' s ' s 



(4) 



where A(a, b, c) = \J 'a 2 + b 2 + c 2 — 2(ab + be + ca). Thus, the masses can be determined by 
measuring these endpoints. Since the distribution is quite steep around the maximum value, 
^UW) snom< i be a good observable to be measured precisely. It becomes 24 GeV even for 
the cases with small mass splitting (I, III), and 80 GeV for the case (II) at ^/s = 250 GeV. 
On the other hand, Ejljjufi measurement may be difficult since the distribution is gradual 
around the minimum, and E^-^ is too small for the cases (I, III) [E^^ = 2.4 GeV at 
v^ = 250 GeV]. 

In the bottom left (right) panel in Fig. [I], we show the Enuj\ distributions for the pa- 
rameter set (IV) with y/s = 500 GeV and £ int = 500 fb _1 . To reduce SM background 



1 The acoplanarity angle </> a cop is denned as the supplement of the difference of azimuthal angles of the 
leptons (jets). 
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FIG. 1. Distributions of dilepton (dijet) energy in dilepton (dijet) plus missing energy event at the 
ILC with yfe = 250 GeV and £ int = 250 fb _1 (top), and with ^ = 500 GeV and £ int = 500 fb" 1 
(bottom). 



contributions, kinematical cuts are applied; aco i > 80° where the acollinearity angle # aco i 
is defined as the supplement of the opening angle of the dilepton (dijet); \M u ^j) — m z \ < 
10 GeV; p^' < ]_gg GeV; |cos#«| < 0.6 for the dilepton case, or |cos0j_,-| < 0.4 and 
|M rec — rriz\ > 10 GeV for the dijet case, where M rec is the recoil mass defined as the invari- 
ant mass of the missing 4-momenta. In the case with the on-shell Z-boson, information of 
the masses can be obtained from the endpoints of the E^{jj) distribution as 

7-,max/min A , a /r\ 

E u(jj) = 1aE ± PaJaP, (5) 
where E = {m 2 A — m 2 H + m|)/(2m^), p = m A /2 • A(l, m 2 H /m\, m 2 z /m\), 7^4 = (s — m 2 H + 
m A )/(2y/srriA) and (3a1a = \fs/{2mA) • A(l, m 2 H /s, m 2 A /s). In spite of the finite-width effect 
of the Z-boson, the endpoints would be seen as sharp edges of the signal plateau. For the case 
(IV), the maximum and minimum of the dilepton (dijet) energies are EJ^-^ = 134 GeV and 
Eyfijft = 181 GeV at y/s = 500 GeV. The signal and background distributions are similar for 

8 



both the dilepton and dijet cases, while the expected number of events are 10 times larger 
for the dijet case than that for the dilepton case. Thus, the statistical errors are small in 
the hadronic signature. On the other hand, the systematical errors would be negligible in 
the dilepton case due to the fine resolution of the lepton momentum measurement. 

For the case with the off-shell Z-boson, from the maximum of the dilepton (dijet) 
invariant-mass distribution, the difference of the two scalar masses can be determined as 

MaSi) =rn A - m H . (6) 

We note that, for the case with small mass splitting, the measurement would be affected by 
QED radiation backgrounds and acceptance cuts. 



B. e + e H + H process 

We here turn to the H + H~ pair production, where H ± predominantly decays into HW ± , 
and W ± further into i^v or qq'. We study the semi-leptonic and all-hadronic decay modes 
as successful signatures. 

First, we study the semi-leptonic decay mode, where the signature is a charged lepton 
plus dijet plus large missing energy. The expected leading background process is r^vjj 
production followed by the leptonic decay of r. The t^vjj background process can be 
reduced by requiring a large recoil mass in the event. The contribution from production of 
fi + fi~jj and missing particles, where one of the muons goes out of the acceptance region, 
also turns out to be negligible. The event simulation for the case (I) [the case (III)] is the 
same as that for the case (II) [the case (IV)], but the overall normalization is multiplied by 
about a half, which is the square of the H ± — > W ± H branching ratio, 68% [73%]. 

In the left and middle panels in Fig. |], distributions of E^ad an d M^d in the semi-leptonic 
decay mode are plotted by using the parameter set (II) at the ILC with y/s = 250 GeV and 
£int = 250 fb _1 with a cut of M rec > 180 GeV. The background contribution is negligible. 

For the case with the off-shell IV-boson, the endpoints of the all-jets (hadrons) energy 
distribution are given by 

(7) 

Here, we note that the invariant mass of all hadrons vanishes at the endpoints. Therefore, 
the hadronic system would be actually observed as one jet near the endpoints. When we 
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FIG. 2. Distributions of E^ad) -^had i n the semi-leptonic decay mode at yfs = 250 GeV with 
Ant = 250 fb" 1 (left and middle) and that of E had at yfa = 500 GeV with Ant = 500 fb" 1 (right). 



apply a cut on the smallest of the dijet invariant-mass at M cut , the endpoints of the energy 
distribution would be replaced by 

nmax/min 771 i /o\ 

^had = lH±E had ± ~1h±PH± Phad, (8) 

with lH± = y/s/(2m H ±), Ph± = (1 - im^/s) 1 / 2 , E had = {m 2 H± - m\ + M 2 ut ) / (2m H ±) 
and phad = rriH±/2 ■ A(l, m 2 H /m 2 H ±, M 2 nt /m 2 H± ). Thus, the mass information can be still 
obtained. Furthermore, the maximum value of the invariant mass of all hadrons is just the 
difference between m#± and m#, 

M™ x = m H ±-m H . (9) 

In the right panel in Fig. [| the -Ehad distribution in the semi-leptonic decay modes 
are plotted by using the parameter sets (II) and (IV) at the ILC with y/s = 500 GeV and 
Ant — 500 fb -1 with a cut of M rec > 150 GeV. Notice that the parameter set (II) corresponds 
to the case where H ± decays into off-shell W and H, and (IV) corresponds to the case where 
H ± decays into on-shell W and H. When the H^-boson is on-shell, the signal distribution 
is like a rectangle where the edges are given by E™^ mm in Eq. (H), but with M cut being 
replaced by mw- 

We note that the dijet system in the semi-leptonic decay mode may be replaced by 
the dijet subsystem which satisfies Mjj ~ in the all-hadronic decay mode where the 
signature is four jets plus missing energy (see below). 
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FIG. 3. Distributions of M rec in the all-hadronic decay mode at yfs = 250 GeV with £; nt = 
250 fb -1 (left), and M rec and M v ; s distributions in the all-hadronic mode at yfs = 500 GeV with 
Ant = 500 fb -1 (middle and right). 



Now, we turn to the all-hadronic decay mode, which results into the signature of four jets 
plus large missing energy. Major SM background comes from the production of four partons 
and two neutrinos. In the left panel of Fig. |3], M rec distribution is plotted for the signal 
process using the parameter set (II) at \fs = 250 GeV with £ int = 250 fb" 1 . To reduce the 
SM background, kinematical cuts of p™ lss > 70 GeV, | cos^m^ < 0.7 and E vis < 120 GeV 
are applied, where # m i SS is the polar angle of the missing 3-momenta and E vis is the sum 
of the energy of all hadrons in one event. As a result, the SM background is sufficiently 
reduced. The minimum of the M rec distribution is at the twice of m#, 

M™ n = 2m H . (10) 

In the middle panel of Fig. ^, the same distributions are plotted but for the signal pro- 
cesses using parameter sets (II) and (IV) at y/s = 500 GeV with £ int = 500 fb -1 . By the 
kinematical cut of Icos^;^! < 0.8, the SM background is sufficiently reduced except at 
M rec ~ rriz- As it is shown in Appendix, the peak of the signal distribution is given by 

M£? = (11) 
m H ± 

It is the advantage of this observable that this relation holds even when the VF-boson in 
H^ 1 — > W ± H is off-shell. Thus, the ratio of m# and can be determined. 

In the right panel of Fig. [3|, the M vis distributions are plotted for the signal processes 
using parameter sets (II) and (IV) at y/s = 500 GeV with £ int = 500 fb -1 . In addition to 
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the kinematical cut applied in the previous panel, the cut of M rec > 110 GeV is applied to 
reduce the SM background which has a missing energy from Z — > vv. After these cuts, the 
SM background is sufficiently reduced except at M vis ~ ra Z - The signal distribution has a 
peak at 

M peak _ mwV^ (m 
1Kis ~ ra H± ' {U) 

when the PF-boson in H ± — > W ± H is on-shell [the case (IV)]. When the W^-boson is off-shell, 
the relation on the peak position no more holds. 



C. Mass Determination 



Here, we summarize the observables for determining the masses of inert scalars. First, we 
consider the determination of m#± and ran in the process e + e~ — > H + H~ . If ra H ± — ran < 
ra w , ra H ± and ra H can be determined simultaneously by measuring the four quantities; 
E^, E^ in Eq. (g), M^f in Eq. (g) and M™ n in Eq. (|TD|). In the left panel of Fig. f§ 
we show how the masses are determined by the measurements of the four quantities for the 
cases (I) and (II) at \fs = 250 GeV. The four bands are plotted on the m#±-m# plane by 
assuming that the four quantities are measured in ±2 GeV accuracy (without any systematic 
shifts). For this assumption, the accuracy of the m#± (ran) determination would be ±2 GeV 
(±1 GeV). On the other hand, if m H ± — ran > mw, the four observables, E^', E^ in 
Eq. (|) but with replacing M cut by m w , M r p e e c ak in Eq. fllj) and M^f in Eq. (pi) are utilized 
for the mass determination. In the right panel of Fig. [|, the four bands are plotted on the 
mn±-m,n plane by assuming that the four observables are measured in ±2 GeV accuracy. 
It turns out that the constraints from measurements of M, p ^ ak and M£f k are more stringent 

Via IcL ^ — ' 

than those from the E™^ mm measurements, if these quantities are measured in an equal 
accuracy. It is expected that peak positions can be precisely determined more than endpoints 
of distributions in the presence of the resolution of energy measurements and the remaining 
background contributions. By combining the four measurements with the uncertainty of 
±2 GeV, mn± and ran can be determined in ±1 GeV accuracy. 

Next, the determination of ra^ can be achieved by combining the observables in the 
process e + e _ — > HA. For the cases with the off-shell Z-boson (I, II, III), EJ^-^ in Eq. (|]) 
and M™*^ in Eq. (|J) measurements can be utilized. However, at -y/s = 250 GeV and 
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FIG. 4. Determinations of rn#± and ran by the four observables are illustrated in the left [right] 
panel for the cases (I, II) [(III, IV)] at yfs = 250 GeV [500 GeV]. Each observable is assumed to 
be measured in ±2 GeV accuracy. 



\fs = 500 GeV, since the two constraints are very similar, these masses cannot be determined 
at one point. In that case, one needs the value of fixed in the process e + e~ — > H + H~ as 
an input to determine m^. While for the case with the on-shell Z-boson (IV), measurements 
of E™*(™ m in Eq. @ can be utilized. In that case, the expected accuracy of the mass 
determination is ±3 GeV for the measurement of the observables in ±2 GeV accuracy. 



IV. DISCUSSIONS 

At the ILC, another important measurement which we have not discussed is the spin 
measurement of the extra particles. In contrast to H and (H, A) in the IDM, supersym- 
metric model contains charginos and neutralinos which are fermions, and the littlest Higgs 
model with T-parity or models with extra-dimension scenario predict heavy M^-bosons and 
heavy photons which are vector particles. Since the production and decay of these parti- 
cles can mimic the signatures studied in this letter, it is particularly crucial to discriminate 
models by the spin measurement of new particles. The spin of new particles could be traced 
in the threshold behavior of the production cross sections, the scattering angular depen- 
dence, and the decay angular dependence. It has been discussed that these measurements 



are possible for the processses we considered at the ILC [^7[ |36|, [37|. On the contrary, those 
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measurements at the LHC are difficult, because the center-of-mass energy in the partonic 
scattering is not fixed and also the center-of-mass system of the partonic scattering cannot 
be reconstructed. 



Finally, we comment on the total comparison for the IDM studies at the LHC and the 
ILC It is shown that the LHC has discovery potential for the inert scalars by using the 
signatures with multileptons plus large missing transverse momentum ||, ^T], |22|, if the mass 
difference is sufficiently large. Thus, in some preferable situations, the evidence of the inert 
scalars would show up at the LHC before the ILC experiment starts. On the other hand, 
as we have shown, the ILC can probe the case with smaller mass difference still in a good 
accuracy, as long as the production processes are kinematically accessible, i.e. y/s > 2m H ± 
and mjj + m-A- Thus, there is a parameter region where any evidence can not be observed 
at the LHC, but the ILC could find it. In any case, at the ILC, the masses of the inert 
scalars can be reconstructed in good accuracy, and furthermore, the spin of the particles can 
be measured in principal. Therefore, the final discrimination of the model beyond the SM 
would be performed at the ILC. 



V. CONCLUSIONS 



To conclude, we have studied the collider phenomenology of the inert scalar particles at 
the ILC in the IDM. The model contains a scalar dark matter candidate which is stable 
due to the unbroken Z2 symmetry, in addition to another neutral scalar and charged scalar 
bosons. At collider experiments, because the dark matter would escape from the detector, 
the signatures always include large missing energy. We have studied the collider signatures of 
the pair production of the charged scalars and the associated production of the two neutral 
scalars with appropriate kinematical cuts to reduce the SM backgrounds. We have also 
investigated the observables to determine the masses of the scalars in these processes. We 
have shown that distinctive signatures can be observed even if the cases with small mass 
difference by applying simple kinematical cuts, and that by combining these observables the 
inert scalar masses can be reconstructed in a good accuracy at the ILC. 
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APPENDIX: INVARIANT-MASS DISTRIBUTIONS INI-> AA ->■ BCBC 

Here, we consider the kinematics of a general process X — > A X A 2 — > BxC\B 2 C 2 , where 
X is a certain initial-state with the fixed collision energy, A { are scalars with the same 
mass rriA, e.g. the identical particles or the charge-conjugate particles. Bi and Ci are any 
particles with their masses mg and mc, respectively, which are produced from the isotropic 
decay of Ai. For the meantime, we consider the case where all the particles are on-shell, i.e. 
y/s > 2m A > 2m B + 2mc- We consider the invariant-mass distribution of the B\B 2 pair: 



where d<& n is the n-body phase-space volume element and Ta is the total decay width of A. 
Using the narrow width approximation for A, it is calculated as 



where (9i, 0) and 9 2 are the decay angles in the rest-frame of Ai and A 2 , respectively. The 
azimuthal angle in the decay of A 2 is fixed to be zero, and the scattering angles in the 
process X — > AiA 2 , which are irrelevant to Mbb, are integrated out. By using the above 
integration variables, (ps 1 + Pb 2 ) 2 is expressed as 




(AI) 



[(p Bl + m) 2 - ^A? + ™ a T 2 a ' l(PB 2 + Pc 2 ) 2 ~ m\Y + m\Y\ ' 




(A2) 



(Pbi +Pb 2 ) 2 



2m 2 B + ![(! + 0a) ( e2 + + 2(3 A eP( Cl + c 2 ) + (1 - P 2 A )P 2 s 1 s 2 c <p ] . (A3) 
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Here, s = p 2 x , /3 A = a/1 - 4m A /s, e = 1 + m 2 B /m 2 A - m 2 c /m 2 A , = X(l,m 2 B /m%m 2 c /rn A ) 
with X(a,b,c) = \J a 2 + b 2 + c 2 — 2(ab + be + ca), Ci = cos^ for i — 1, 2 and = cos</>. 
Although the numerical integration in Eq. ( [A2] ) is straightforward, we find that, unless 
y/s ~ 2771^4, it is a good approximation to neglect the last term in Eq. ( |A3| ) . In that case, 
the Mbb distribution can be analytically expressed as 



da 



dM 2 BB 



oc log 



A{(s~2m 2 A )(M 2 BB -2m 2 B )-m\e 2 } 



for M min < M BB < m By /s /m Al and 
da 



dM 2 BB 



oc log 



s i{(l + p 2 A )p-2(3 A e} 
s 2 {(l + p 2 A )p + 2p A e} 



(A4a) 



(A4b) 



4{(s-2m 2 A )(M 2 BB -2m 2 B )-m\e 2 }_ 
for mB\fs/tTiA < Mbb < ^fmax, where M max / min = \/s/2- (e±(3 A (3). Thus, we find that the 
distribution has a peak at mB\/s/mA, which therefore depends only on the ratio ms/mA 
but not on mc- 

Here we give some comments. In case with rriA < itib + mc, if the particles Cj are off- 
shell, the peak position of the M BB distribution is still given by \/sm B /mA- On the other 
hand, if the particles Bi are off-shell, the reconstructed M BB distribution is largely deformed. 
The Mac distribution can be obtained by replacing B -H- C from the Mbb distribution in 
Eqs. (|A4[). Thus, the peak position of the Mcc distribution is given by M^ k = mc\/s/mA, 
independently from tub- 

At e + e~ colliders where the total 4-momenta of the collision can be fixed, it is possible to 
assess the invariant mass of the all missing particles (the recoil mass) by M 2 ec = (p in — p v is) 2 , 
where p in is the total 4-momenta of the initial e + e~ system. Therefore, the invariant mass of 
the pair of two missing particles can be reconstructed, if these are the only missing particles 
in the event. 
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